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Abstract 

We show that the physical minimum of the Constrained MSSM is only free from dan- 
gerous charge and colour breaking minima in the region of parameter space bounded by 
110 ;$ fnx/2 ~ 400 GeV and 80 mo < 170 GeV. In the remaining regions the cosmology 
is severely constrained. 
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1 Introduction 



Unphysical charge and colour breaking (CCB) vacua (T|, ||, ||, £|, §, §, [7], § have come under 
renewed scrutiny recently^. For a number of models it has been found that CCB vacua are 
present in the whole of the parameter space which has not already been excluded by exper- 
iment. This is true for models where supersymmetry breaking is driven by the dilaton [[|, 
for M-theory in which supersymmetry breaking is driven by bulk moduli fields |7], |8| and for 
the MSSM at the low tan/5 fixed point f|. 

In this letter we extend the analysis of Ref. || to a complete determination of CCB 
bounds in the Constrained MSSM (CMSSM). The present work is partly an update of the 
results of Ref.||, and from that study and also Ref.|| it is clear that the parameter space 
is bounded from different directions by CCB bounds, dark matter bounds and experimental 
bounds. Dark matter bounds (the requirement that neutralino dark matter be consistent 
with a universe that is 12 Gyr old) tend to eliminate regions where mo is large. On the 
other hand, experiment (in particular Higgs and chargino searches) is eliminating regions of 
low rrii/2- CCB bounds however tend to 'favour' regions where mi/ 2 ^m or in other words 
high m and low m^j. What we shall demonstrate in this letter is that there is only a 
small amount of remaining parameter space in the CMSSM which does not have global CCB 
minima. 

Before continuing, we emphasise that there are various schools of thought regarding CCB 
minima. This is because the tunneling rate from the physical vacuum into any global CCB 
minima is extremely small, so that the physical vacuum is essentially stable on the lifetime 
of the observable universe. The authors of Ref. 0, for example, invoke a principle that 
'the cosmological constant is zero in the global minimum' in order to explain the vanishing 
cosmological constant, thus requiring that the standard model minimum be deeper than any 
CCB minima. If this were the case, then CCB minima would indeed be a severe problem 
and would impose severe constraints on the MSSM. By contrast, some other authors take 
a minimalist approach; if there is no chance of tunneling, then CCB minima are not a 
problem. However, one has to explain how the physical vacuum was chosen over the much 
wider CCB minima [|], ||. In fact this does naturally occur in some models, but not in 
others (e.g. supergravity models which possess a Heisenberg symmetry O, including no- 
scale models of supergravity [Kj). As the vacuum choice depends on unknown details of 
our cosmological history (e.g. the inflationary potential) we think that CCB minima should 
ultimately be regarded as a constraint on early cosmology rather than particle physics. We 
also emphasize that since tunneling between vacua is so slow as to be irrelevant, it is more 
appropriate to flag models which contain any CCB minima, regardless of whether they are 
global or local. 

We begin by describing the most dangerous directions, assuming the usual i?-parity 
invariant superpotential of the MSSM, 

Wmssm = h v QH 2 U c + h D QHiD c + h E LH x E c + ^H^, (1) 

1 In Refs.|| |tJ the bounds which will be of most interest here were, in deference to historical precedent, 
referred to as Unbounded From Below (UFB) bounds. This is a confusing misnomer since the directions are 
not, in most cases, unbounded from below, and we revert to calling them CCB bounds. 
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and a degenerate pattern of supersymmetry breaking with universal scalar mass 2 parameters 
(wig), trilinear couplings (A ) and gaugino masses (m^) at the GUT or Planck scale. The 
notation is the same as in Ref. ||. 

As emphasised in Ref. 0, the flat directions in the MSSM are a direct result of the 
adoption of i?-parity, to prevent the proton decaying. The dangerous F and D flat directions 



which will be of interest in this letter are constructed from gauge invariants involving H 2 \11 



12] . This is because its mass squared parameter, m 2 ., appears in the potential along these 



directions, and it must be negative in order to drive electroweak symmetry breaking. The 
first example of this kind in the literature is (see Komatsu in Ref.[|l|]) 

LiQ 3 D 3 ; H 2 Li (2) 

where the suffices on matter superfields are generation indices. With the following choice of 
VEVs; 

h° 2 = -a 2 /j,/h D:)3 
di 3 = dji 3 = ajj,/ho 33 

i>i = a\/l + a 2 n/hr) 33 , (3) 

the potential along this direction is F and .D-flat, and depends only on the soft supersym- 
metry breaking terms; 

V = _^ a 2 (a 2 (m 2 + m 2j + ^ + ^ + m 2 Q J_ (4) 

At large values of a 3> 1 the potential is governed by the first term. Because m 2 is required to 
turn negative during the renormalisation group running (for successful electroweak symmetry 
breaking) the potential can develop a charge and colour breaking minimum at a scale of 
few x iijhvy Ensuring that this does not happen leads to the constraint in which we are 
interested. 

The above is not quite, but is very close to, the deepest 'fully optimised' direction [|]. H- 
To get the fully optimized condition we parameterize the VEV of U{ as [I | 

h = ^LCi 2 ^/h D33 . (5) 

Minimisation of V with respect to 7l then gives 



2 1 + a 2 2m 2 ^ 2 
^ = — — - 1^±2~ ( 6 ) 



2 _ l -r u _ ^" b Li 

a 2 g 2 a A fi 

j2 



where g = (g' + g%)/2, and a potential of 

2 4 
U D 33 9 

When 7 2 < then one should set 7! = to get the most stringent condition although, for the 
regions of parameter space of interest here, this will never be the case. The potential typically 
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has a depth of >, 10 6 m^ at the minimum. In addition, as shown in Ref. 0, the depth of 
the minimum is actually extremely sensitive to the choice of soft supersymmetry breaking 
parameters, so that, for most reasonable choices of parameters, the difference between Eq.(|4]) 
and Eq.([7D is negligible. In Ref.|§, emphasis was also placed on the closeness of the usual 
condition (that the physical vacuum be the global minimum) to the more relevant condition 
(that the physical vacuum be the only minimum). There is generally an extremely thin, but 
cosmologically interesting, region of parameter space between the 'allowed' and 'disallowed' 
regions, where there is a CCB minimum which is local. 

In order to obtain the bound we now need to take account of the renormalisation group 
running of the mass-squared parameters between the weak and GUT scales. To do this we 
shall assume that the largest mass, and therefore the appropriate scale at which to evaluate 
the parameters is <fi = hu Z3 (h!l). This minimises the top quark contributions to the effective 
potential at one-loop. Further corrections to the potential are assumed to be small. As shown 
Ref. 0, this approximation is adequate for determining CCB bounds on the supersymmetry 
breaking parameters (although they also note a curious bifurcation in behaviour when one 
varies the scale at which the parameters are evaluated). 

In the above potentials, (h®) = —a 2 ^/h£ )s3 ^E ss so that the Eq.(|4]) is of the form 

V = ^MA + B/b) (8) 

where A = ra^ifi) + m| .(0), B is the other combination of mass-squared parameters (also 
evaluated at 0) which appears in the potentials above, 

= 0/ilW (9) 

and 

m = M GUT h D33 (1Q) 
h U33 fi 

for the LQD, LH 2 direction described above, or 

m = M GUT h E33 (n) 
h U33 fJ, 

for the equally dangerous LLE, LH 2 direction. The traditional (no global CCB minima) 
bound is saturated by V — V — 0; the non-trivial solution is therefore also a solution to 
V — V — where 

V = (f)A + B/b. (12) 

Hence, only the GUT scale parameters plus the parameter b enter the traditional bound as 
was pointed out in Ref. || . In fact the bound always becomes more restrictive with increase 
in b, since this decreases the positive contribution to the potential. 
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2 The analysis 



2.1 CCB Minima 

First we shall examine the properties of the CCB minima when we vary the soft super- 
symmetry breaking parameters. We will not however calculate the tunneling rates which is 
generally a complicated numerical task. In this instance, it is made particularly difficult by 
the presence of more than one gauge invariant involved in the flat directions. (As an aside 
we note that, for the directions involving only the single UDD operator examined in Ref. ||, 
a very accurate analytic approximation can be found for the tunneling rate.) However, from 
dimensional considerations and estimates of tunneling rates in Refs.[§], §], it is clear that the 
lifetime for these directions will be much longer than the age of the universe. 

We can make additional observations on the stability of physical vacuum by adopting 
the approach of Ref. || in the light of Ref. . In the latter it was pointed out that the CCB 
bounds are least restrictive when A = —7711/2 Q- Hence in Fig. [l| we show contours of the 
VEV of the CCB minimum in the mo, 7774/2 parameter space, for tan (3 = 2, 10 and fi < 0, and 
taking A = — 7711/2 and 0. Here we have plotted log 10 (VEV/246 GeV) for the LLE, LH 2 
direction. Above the solid contour there are no global CCB minima. In the thin shaded 
region (approximately 10 GeV wide) above the solid contour, there are CCB minima which 
are only local, and the physical minimum is still the global minimum. In Ref.0 a general 
survey was made of all the possible dangerous directions in field space, and this direction 
was found to give the severest bounds (at least for the CMSSM) so henceforth we shall 
only be considering this. The figure illustrates that the VEVs of CCB minima are smallest 
when the CCB bound is close to being saturated. This leads to the rather counter-intuitive 
fact that the physical vacuum is least stable in these regions. When the bound is strongly 
violated the VEVs become very large and the lifetime of the physical vacuum increases. 
The CCB bounds are least restrictive when A « —7771/2, and they are close to the analytic 
approximation found in Ref.0; if we define mo = mo/m^ then the analytic approximation 
is given by 

~ 2 . /(3m§) -fl(3rrag)(l - p p ) 

m ° > I^3i P [l6) 

where 

f(x) = 1.43 - 0.16x + 0.02x 2 
g(x) = 2.94- 0.2x + 0.02x 2 

l/p p = l + 3.17(sin 2 /3-sin 2 /3 QFP ). (14) 

These functions were evaluated in Ref. 0] for mi/2 = 200 GeV and in the one loop approxi- 
mation. However we use the full two loop value of ta,n(3® FP which is ~ 1.6. (Note that the 
important factor here is the distance from the fixed point which is given by p p . For a given 
value of tan (3 this obviously depends sensitively on the fixed point value, tan j3 QFP , so we 

2 In our sign conventions, the 3-4 element of the neutralino mass matrix is — fi, and the mixing term in 
the stop mass matrix is m t [A t + fi cot j3) . 
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cannot expect the analytic approximation to be better than ~ 15%.) For tan/3 = 2, 10 and 
A = —mi/2 we find the bounds m ^,130,40 respectively. The first compares favourably 
with the bounds at m-1/2 = 200 GeV found numerically in Fig. |2| but the tan /3 — 10 bound is 
found to be larger than the analytic estimate. This suggests that the bottom quark Yukawa 
and/or the two loop contributions to the beta functions are already contributing significantly 
to the bounds at tan/3 = 10. 



2.2 Experimental Constraints 

We now discuss the experimental and cosmological bounds on the CMSSM parameter space. 
Recent runs at LEP at center-of-mass energies of 172 and 183 GeV have excluded large 
areas of the CMSSM parameters space, and subsequent runs at ~ 190 and 200 GeV will 
push the bounds even further. In the CMSSM, the dominant constraints at moderate to 
high tan/3 come from searches for chargino pair production, and modulo a small loophole 
which can occur when the mass of the sneutrino is close to the chargino mass, the experi- 
mental bounds saturate the kinematic limit of m x ± ~ 91 GeV . Chargino iso-mass contours 
of 91 and 100 GeV , representing the current and projected LEP 200 chargino mass bounds, 
respectively, are displayed in Fig. |2| as dashed lines in the m.1/2 — m plane, for two represen- 
tative values of tan /3 and both signs of //. Also shown as a dotted line is the current LEP183 



slepton mass bound [13|], which is roughly 84 GeV for large mj — m^. 

At low tan /3, the dominant CMSSM constraint comes from searches for Higgs production 
at LEP183. Not only are the experimental bounds strongest at low tan/3, roughly 90 GeV 
at tan/3 = 2, but here the tree level Higgs mass is also smallest (m^ ree m rriz cos 2/3 for 



rriA 3> rnz)- Radiative corrections to the Higgs mass fl4 |, which depend logarithmically on 
the sfermion masses, must then be very large, leading to strong lower bounds on the masses 
of the sfermions, and in particular the stops. However, the extraction of the radiatively 
corrected Higgs mass in the MSSM has an uncertainty of ~ 2 GeV, so we conservatively take 
rrih > 88 GeV as our experimental lower limit at low tan/3. The LEP183 Higgs bounds are 
shown in Fig. ||] as a dot-dashed line. The entire displayed region of Fig. ^|a is excluded by 
the Higgs mass constraint, while none of the displayed region in Fig. ^|d is in conflict with 
the current Higgs bound. Of course the Higgs mass corrections are very sensitive to tan /3 for 
tan/3 near the quasi-fixed point, and, additionally, the experimental limit falls for tan/3 > 2, 
and so the Higgs constraint moves quickly to the left for tan/3 > 2. The chargino bound, 
on the other hand, moves to the right, and for 100 GeV < m < 200 GeV, the two bounds 
together exclude 7711/2 £ 110 GeV for all tan/3 and both signs of //. 



2.3 Cosmological Constraints 

Over most of the CMSSM parameter space, the lightest supersymmetric particle (LSP) is 
a bino-like lightest neutralino \. .R-parity ensures that the LSP is stable over cosmological 
time scales, and in the CMSSM, the LSP typically has a cosmologically interesting relic 
density [|l^, [16[]. In much of the parameter space, in fact, the relic abundance of neutralinos 
is so large that it is in conflict with the observed age of the universe, tjj > 12 Gyr, and 
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the corresponding upper limit of Q x h 2 < 0.3 can be used to exclude large areas of m and 
mi/2, as follows. In the early universe, gaugino-like neutralinos annihilate predominantly 
via sfermion exchange. Increasing either ttlq or mi/2 drives up the sfermion masses, lowers 
the annihilation rate, and increases the neutralino relic abundance. Thus an upper bound 
on Q x h 2 translates into an upper bound on the parameters ttt-o, 777,1/2 ■ The experimental and 
cosmological constraints are nicely complementary, in part because the cosmological bounds 
put upper limits on the parameters for which particle searches provide lower limits. 

The dark shaded area in Fig. ^ delimits the cosmologically preferred region with 0.1 < 
Q x h 2 < 0.3. The upper limit, as described above, comes from an upper bound on the age 
of the universe. The lower limit is more of a preference than a bound, stemming from the 
wish to have the neutralinos comprise a significant fraction of the dark matter. The two 
narrow vertical channels arise from s-channel neutralino annihilation on the Higgs and Z° 
poles (for some values of tan/5 they have merged into one large pole region). Note that the 
top of the shaded area intersects with the "theory" excluded area (where a stau is the LSP) 
at 777,1/2 ~ 450 GeV . Thus £l x h 2 < 0.3 yields an upper bound on 777,1/2 and on m (except in 
the close vicinity of the pole region, which is largely excluded by the LEP chargino searches). 
At sufficiently low tan /3, the Higgs bound moves to the right of the shaded region, and the 
incompatibility of a Higgs lower bound of 86 GeV with the lower limit on the age of the 



universe excludes values of tan/3 less than 2, for \i < 0, or 1.65 for /i > [[17] . 

The dark solid line in Fig. [I] is reproduced in Fig. |^, where we have chosen Aq = —7771/2 to 
minimise the size of area containing CCB minima (see above and Ref. M). We note that there 
is only a restricted region in {7711/2, mo} which is both cosmologically and experimentally vi- 
able and which is free of CCB minima, and moreover, such a region exists only for tan (3 > 2.3 
for fx < and tan/3 > 2.1 for /j > 0. At large tan/3 >, 20, s-channel annihilation through the 
pseudoscalar Higgs and heavy Higgs can contribute significantly to the neutralino annihi- 



lation cross-section, and the above cosmological bounds are weakened ||18|| . Coannihilations 
with light sleptons also substantially reduces the relic abundance of neutralinos near the line 
rrif = m x and can provide a window at large 7721/2 [[19]. 

3 Conclusions 

In this paper we have re-examined the unphysical charge and colour breaking minima in 
the Constrained MSSM. In summary, we find that most of the parameter space which is 
not already excluded by experiment or by cosmological considerations has an unphysical 
vacuum which is lower than the physical one. The region of parameter space which remains 
is roughly 

110 < mi/2 < 400 GeV 

80 < m < 170 GeV, (15) 

although there is a narrow region at 7711/2 ~ 150 (110) GeV for tan/3 = 10(3) for which 
higher values of m are still allowed by the projected chargino searches. Future tri-lepton 
searches at the Tevatron should push the lower bound on 7771/2 to from 180 — 240 GeV for 



the regions of Fig. |] |20| and will help close this loophole to larger mo- 
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The regions which do have an unphysical minimum, as we stated in the introduction, are 
not completely excluded but they must have a constrained cosmology. Some general points 
were made in Refs. 0, so let us apply some of these considerations to this specific model. 
Our main requirement for any acceptable cosmological scenario is that it should explain why 
the physical minimum is chosen instead of the global minimum. There are currently three 
possible explanations; 

• A high reheat temperature after inflation 

• Extra contributions to the effective potential during inflation 

• Breaking of i?-parity 

In the case that the field is driven to the origin during reheating, there is an additional 
constraint in all models in which supersymmetry breaking is transmitted to the physical 
sector gravitationally. In these models the gravitino mass is typically of the same order as 
the weak scale, and successful nucleosynthesis requires that 

T reheat < 10 9 GeV. (16) 

In order to lift the unphysical minimum we should have a reheat temperature which is greater 
than roughly 10 _1±1 x the scale of the VEV (the 'uncertainty' coming from the shape of the 
effective potential and the contribution of light particles to it). Hence a particular T re h eat can 
be effective in the region bounded by the corresponding contour in Fig. (|ip. From Fig. (|I|) 
we see that the condition (^) can quite easily be satisfied away from the fixed point in the 



CMSSM. 

If the physical vacuum is chosen because of extra terms generated during inflation then 
there are constraints on the possible inflationary potentials. For example potentials with a 
Heisenberg symmetry are eliminated in this case because they are flat at tree level and the 
one loop corrections to the mass-squareds tend to be negative || [10 . 



If .R-parity is broken explicitly [EL |6| we must choose a symmetry other than .R-parity 
to prevent the proton from decaying, and of course lose the neutralino as a dark matter 
candidate (but, on the plus side, free up a large region of parameter space in Fig. (H)). 
Note that, as well as the lepton number violating version of -R-parity violation considered in 
Ref.@, the CCB minima would also be lifted by a large Majorana mass for a right handed 
neutrino. If i?-parity is broken during a stage of pre-heating then it is in principle 
possible to have neutralino dark matter which is stable. It is also possible that further lifting 
of the potential would occur in any case from extra operators required in the visible sector. 
Both of these mechanisms work because at large field values (i.e. those corresponding to the 
VEV of the unphysical minimum), the MSSM is no longer a good description of the effective 
potential, and the potential is no longer F flat along the relevant directions. 
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Figure 1: The location of the CCB VEV, in the direction LLE, LH 2 , for tan/5 = 2, 10 and 
A = 0, — mi/2- The dashed lines are contours of constant log 10 (VEV/246 GeV). In the 
shaded strip, the CCB minimum is not global. There are no CCB minima above the shaded 
region. 
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Figure 2: The combined cosmological and experimental constraints on the constrained 
MSSM, for tan (3 = 2, 10 and both fi < and \x > 0. The dashed contours represent current 
and future LEP chargino bounds, dotted contours are slepton bounds, and dot-dashed con- 
tours are Higgs bounds. The light-shaded region gives 0.1 < £l x h 2 < 0.3. Below the solid 
contour, CCB minima are present in the LLE, LH 2 direction. We have chosen A = -m 1( / 2 
to minimise the area containing CCB minima. 
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